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Mean encounter times for cell adhesion in hydro dynamic flow: 
analytical progress by dimensional reduction 
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Abstract. - For a cell moving in hydrodynamic flow above a wall, translational and rotational 
degrees of freedom are coupled by the Stokes equation. In addition, there is a close coupling of 
convection and diffusion due to the position-dependent mobility. These couplings render calcu- 
lation of the mean encounter time between cell surface receptors and ligands on the substrate 
very difficult. Here we show for a two-dimensional model system how analytical progress can be 
achieved by treating motion in the vertical direction by an effective reaction term in the mean first 
passage time equation for the rotational degree of freedom. The strength of this reaction term can 
either be estimated from equilibrium considerations or used as a fit parameter. Our analytical 
results are confirmed by computer simulations and allow to assess the relative roles of convection 
and diffusion for different scaling regimes of interest. 
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Introduction. — Biological function is often based on 
the formation of a specific binding complex between recep- 
tor and ligand [1]. However, in order for binding to oc- 
cur, a physical transport process must exist which brings 

' the binding partners to sufficiently close proximity [2] . In 
many cases of interest, this transport process is rather 

' complex. Usually it contains several coupled degrees of 
freedom, like a cell surface receptor moving laterally on 
a membrane which fluctuates in the vertical direction [3] . 
The efficiency of biological transport processes often can 
be framed as mean first passage time (MFPT) problems, 
for example for the gating of ion channels [4] or the ar- 
rival of a virus at the nucleus [5]. Another example of a 
complex transport process of large biological relevance is 
the receptor-mediated adhesion of cells which are carried 
over a ligand-coated substrate by hydrodynamic flow [6]. 
Here the mean encounter time between receptors and lig- 
ands is a measure for the efficiency of cell adhesion under 
the conditions of hydrodynamic flow [7]. For this system, 
additional complications arise from the presence of multi- 
ple length scales. For the micron-sized cell, the hydrody- 
namic equations result in coupling of the translational and 
rotational degrees of freedom. Even for high shear rates, 
Brownian motion is relevant because receptors and ligands 
are nanometer-sized objects, thus even small movements 
for the cell result in a large effect on the molecular level. 



Here we show that despite the presence of these complica- 
tions, analytical progress can be achieved by dimensional 
reduction of appropriate degrees of freedom. The reduced 
description then contains effective parameters which have 
to be obtained from the full model. 

Experimentally the binding of cells to a substrate in hy- 
drodynamic flow is often studied in flow chambers because 
this setup allows for controlled flow conditions [8,9]. In 
vivo, this situation is relevant for white blood cells, which 
travel the body with the blood flow, but have to adhere at 
very specific locations, e.g. close to sites of inflammation. 
Similar mechanisms are used by cancer and stem cells. 
Moreover, malaria-infected red blood cells also undergo 
adhesion to the vessel walls under flow conditions. Apart 
from using flow chambers, one can further reduce the ex- 
perimental system by employing biomimctic analogs of 
cells, that is receptor bearing micro-beads [10,11]. The ef- 
ficiency with which cells or beads in flow can bind to a sub- 
strate depends crucially on the spatial distribution of re- 
ceptors and ligands [12]. Previously, we proposed a model 
based on a Langevin equation for a spherical particle in 
linear shear flow above a wall which allows to numerically 
compute MFPTs for different ligand and receptor dis- 
tributions and flow parameters both for two-dimensional 
(2D) and three-dimensional (3D) movements [7,13]. Due 
to the complex geometry arising from the receptor and 
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ligand distributions and the complexity of the position- 
dependent mobility functions arising from the hydrody- 
namic equations, exact analytic results for the MFPT can- 
not be obtained in this general case. 

In this letter we present analytical advances for 2D mo- 
tion with homogeneous ligand coverage on the substrate. 
In our model system, motion can occur only in x- and 
z-directions and rotation is restricted about the y-axis. 
Due to the assumption of homogeneous ligand coverage, 
motion in the x-dircction is not relevant. Thus we deal 
with two degrees of freedom, falling in z-direction and ro- 
tation about the y-axis. This is the simplest model sys- 
tem which combines rotational and translational degrees 
of freedom in a non-trivial manner. Because cell move- 
ment usually occurs in the regime of small Reynolds num- 
ber, their coupling is determined by the Stokes equation 
for viscous flow. In addition, we account for Brownian 
motion which is ubiquitous in biological systems on the 
nanoscale and essential for receptor-ligand binding to oc- 
cur. Here we show that this model system can be further 
reduced by effectively integrating out the translational de- 
gree of freedom. This results in an ordinary differential 
equation for the MFPT of the rotational degree of freedom 
which includes a non-trivial reaction term that represents 
the falling motion of the cell. We show that this equation 
can be solved analytically. From this solution we then 
derive various expressions that describe asymptotic lim- 
its. By comparing with computer simulations, we finally 
show for which parameter range our analytical results are 
valid. Moreover, our analytical calculations still explain 
the numerical results outside this parameter range when 
the effective reaction rate is used as a fit parameter. 

Model definition. — We consider a sphere of radius 
R moving in linear shear flow with shear rate 7 above a 
planar wall. As explained above, we restrict its motion 
to two dimensions, that is the translational motion of the 
sphere is restricted to a plane perpendicular to the bound- 
ary wall, i. e. the xz-plane, and rotations are only allowed 
about the y axis (see fig. la). As depicted in fig. la the 
circumference of the sphere lying in the xz-plane is cov- 
ered with N r equidistantly distributed receptor patches 
of height r <C R and radius r p <C R. The boundary 
wall is homogeneously covered with ligands. In order to 
drive it onto the substrate, in vertical direction the sphere 
is subject to a constant force — F z . In experiments, this 
force arises from gravity because cells or micro-beads are 
usually slightly denser than the surrounding medium. A 
receptor-ligand encounter occurs with certainty whenever 
a receptor patch has some overlap with the boundary wall. 
Because we also consider Brownian motion, the receptor- 
ligand encounter is stochastic. Our goal is to calculate the 
corresponding MFPT. For the calculation of first passage 
times the motion in x direction can be neglected as the 
system is translationally invariant in this direction due 
to homogeneous ligand coverage. Moreover, the regular 
distribution of receptor patches generates a 6 S := 2w/N r 



symmetry for the 9 coordinate (i.e., the angle about the 
y-axis). Therefore, we deal with the situation illustrated 
in fig. lb of a diffusive particle moving in the # s -periodic 
(z9) plane (i.e., a cylindrical surface with z > R) with ab- 
sorbing boundaries 8(z) (solid lines in fig. lb) representing 
the boundaries of the encounter areas. 

If T(z', 9') is the MFPT to reach a point on the absorb- 
ing boundary when started at some point (z' , 9') then for 
practical purposes only (T(z' ,9'))g> , i.e., the MFPT av- 
eraged over all initial orientations, is a relevant quantity 
as it is experimentally very difficulty to prepare a certain 
initial orientation. Concerning the initial height we previ- 
ously showed [13] that for z' > z m > R + r the following 
relation holds true 



{T(z',9')) e ^T(z m \z') + {T( 



(i) 



where T(z m \z') is the mean time to fall from the initial 
height z' to the intermediate height z m and 9 m is the 
sphere's orientation at height z m . Eq. (1) states that 
if the angle averaged MFPT is known for some height 
z m the MFPT can be calculated for any other height 
z > z m without further considering the orientational de- 
gree of freedom. Moreover T(z m \z') can be calculated 
exactly [13]. Thus for the following we choose as the 
initial height of the disk always z' = R + r and write 
T(9') :=T(z' = R + r ,6'). 

Model reduction to one dimension. In order 
to calculate (T(9 / ))gi one must in principle consider the 
motion in the full (z6>)-plane. Here, we show that under 
certain conditions a good approximation formula for this 
MFPT can be derived by considering only the motion of 
a particle in a periodic one-dimensional 0-space. For that 
we first approximate the area in the (z0)-plane in which 
encounter occurs by a rectangle of width 2# (see fig. lb). 
As an appropriate 6*o we choose the mean width of the 
encounter interval the particle sees while being at height 
z < R + ro , that is 
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dz9 (z)^ s (z), (2) 



where 
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-F z (z-R)/k B T a 



(3) 



is the stationary probability distribution for the height of 
the sphere z with Boltzmann's constant ks and ambient 
temperature T a . 9q{z) is the boundary of the encounter 
area and depends on receptor height and radius through 
6*o (z) = &vccos(z/{R + r )) + r p [13]. Thus 9 increases 
with r p . 

Changes in 9 are due to rotational diffusion with dif- 
fusion coefficient Dg and a drift Ag oc 7 arising from the 
linear shear flow. This suggests to take the corresponding 
Fokker-Planck equation as a starting point for a reduced 
model. In order to account for the effect of the motion 
in z-direction, we argue that for 9 G [0, 28 ] an encounter 
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Fig. 1: (a) A sphere (radius R) in linear shear flow (shear rate 
7) covered with iV r = 6 receptor patches with height ro and 
radius r p is shown. The sphere translates in the zz-plane. The 
orientation about the symmetry axis is given by the angle 8. 
(b) The encounter problem can be mapped to that of a particle 
moving in z6 space, where the motion in ^-direction is 9 S := 
2n/N r periodic as the receptors are equidistantly distributed. 
The solid curves shown enclose the areas in which receptor- 
ligand encounter occurs. For the calculation we approximate 
these areas by a rectangular area (grey shaded) of height ro 
and width 20o- 



occurs only with a finite probability. Thus, for the prob- 
ability p(9, t) to have the orientation 9 at time t we setup 
the following reaction-advection-diffusion equation for the 
interval [0, 9 S [ 

d tP (9, t) = D g d 2 eP {9, t) - Agd eP (9, t) - k(9)p(9, t). (4) 

The reaction term is defined by k(9) = k for 9 6 [0, 2#o] 
and k(9) — otherwise. It accounts for the finite proba- 
bility that encounter occurs while the particle is passing 
the interval [0, 29o}. This approximation makes sense if 9q 
is small and if in addition 9q/9 s <C 1. In our context, this 
assumption is justified because we have the separation of 
length scales r p ,r <Si R- Then the time for passing the 
interval [0, 2#o] is short compared to the time between suc- 
cessive receptor-ligand encounter and the correlation be- 
tween the heights z for successive encounters are small, 
thus justifying our mean field treatment for the transla- 
tional degree of freedom. The assumption 9q/9 s <C 1 also 



suggests to take the limit 9q while keeping 29ak =: k 
constant. Then the term —k(9)p(9,t) in eq. (4) becomes 
—kS(9)p(9,t). eq. (4) is a differential equation for the ro- 
tational motion in which the effect of vertical motion has 
been absorbed into the new model parameter k. 

Boundary conditions. — Because of the 9 S - 
periodicity the probability p(9, t) must fulfil periodic 
boundary conditions, i. e. p(9 + 9 s ,t) — p(9,t). Further- 
more, integrating eq. (4) over a full period 9 S we obtain 
for the total loss of probability 

j t J° +(>S d9p(9 1 t) = -kp(0,t), (5) 

i. e. k denotes the rate of absorption at the boundary. In 
the limit k — > 00 we have purely absorbing boundaries 
with p(0, t) = p(9 s ,t) — 0. On the other hand integrating 
eq. (4) over the open interval ]0, 9 S [ we obtain 

d r e '~ e 

— lim / d9p(9, t) = J(0+, t) - J(8 s -,t) (6) 
at e^o J 0+e 

with the probability current J(9,t) := —(Dgdg — 
Ag)p(9,t). Combining eq. (5) and eq. (6) we see that the 
reactive delta-function at 9 = is equivalent to so-called 
radiation boundaries [14] 

J{9 s ,t)-J(0,t)=kp(0,t). (7) 

In the limit of zero encounter probability, i. e. k = 0, the 
probability flux leaving the interval on the right boundary 
is equal to that entering the interval at the left boundary 
and the total probability is conserved. 

Mean first passage time. — Let G(9',t) denote the 
survival probability at time t, i.e. the probability that 
no encounter has occurred until time t under the condi- 
tion that the initial orientation of the particle at time 
t' = was 9'. With p(9,t\9',0) being the conditional 
probability for the particle of having the orientation 9 
at time t when the initial orientation of the particle at 
t' = was 9', the survival probability can be written 
as G{9',t) = fH s &9p{9,t\9' ,0). G(9',t) obeys the adjoint 
equation [15] 

d t G(9', t) = {Dgdl + Agdg, ~ kS(9'))G(9' , t). (8) 

The mean first passage time T{9') follows from the survival 
probability as T{9') — J °° dtG(9' , t). Thus, an equation 
for the MFPT is obtained by integrating eq. (8) over all 
times. With lim^o G(9', t) = and G(6»',0) = 1 the 
MFPT is the solution of [15] 

{Dgdl + Agdg, - kS(9'))T(9') = -1. (9) 

The general solution of eq. (9) is given by 

T{e ' ) = -Tg + % (10) 
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with two integration constants ai,a^ which have to be 
determined in order to match the boundary conditions, i.e. 
periodicity T(0) = T{9 S ). The condition corresponding 
to eq. (7) is dg,T\ e , =0 - d 0l T\ e , =6s = (k/D e )T(0), which 
follows from integrating eq. (9) from — e to e with some 
e > and then taking the limit e — > 0. This gives 



T{9') 



Ag k ' 



(11) 



This equation has to be corrected by the physical expec- 
tation that T(9') = for 9' G [0, 26» ] as we used for the 
initial height always z' = R + r . Therefore, for these 
orientations encounter occurs instantaneously. This prob- 
lem is fixed by replacing 9 S by A9 := 6 S — 29q in eq. (11). 
Then averaging over all initial orientations including some 
for which T is zero we arrive at our central result 



(T(e'))e 



i r 

= T J d6'T(6') 

A e A0 2 COth(4^) -2DgAg M 2 



k9 x 



(12) 



As one would expect {T{9'))gi becomes the smaller the 
larger the encounter probability, i.e. the larger the reac- 
tion rate k. 

Parameter estimates. — In order to apply eq. (12) 
to the situation of a particle moving in the (z#)-plane we 
still have to provide expressions for the reaction rate k 
as well as for the diffusion constant Dg and the drift Ag. 
Regarding the reaction rate, we have to consider the full 
system again. We first note that in the stationary state the 
probability for the sphere to be at a height between z and 
z + dz is given by & s (z)dz following from eq. (3). Thus, 
we get for the probability for encounter while the sphere is 
oriented such that 9 G [0, 29 ] (valid in the limit 9 <C 1) 
p z := 1 — cxp(— F z ro/fcsT a ). On the other hand using 
the originally introduced rate k the mean probability for 
encounter is p^ := 1 — exp(— fcr) where r is the mean time 
it takes the particle to pass the interval of length 29q. For 
purely diffusive motion Ag = we estimate r = 9^/2Dg 
which is the mean first passage time to reach the bound- 
ary of [0, 2#o ] when initially started at 9(t = 0) = 6q. In 
the limit of very large drift motion becomes purely deter- 
ministic and we expect r — 29 /Ag. These two limiting 
cases may be combined to provide 1/r = 2Dg/9^+Ag/29o- 
Under the assumption that the position between two suc- 
cessive approaches of the interval [0, 29 ] is independent 
one can get an estimate for k from the condition pj. = p z 



k B T a r k B T a \ 9 +Ae ) 



(13) 



where we used k = 29ok. From eq. (13) one sees that en- 
counter becomes the more probable the larger the vertical 
drift force F z and the faster rotation given by Dg and Ag. 



For a sphere in linear shear flow above a wall the coef- 
ficients Dg and Ag depend on the height z as 



De(z) 



^0 rr (R/z), A (z) = l(l-p dr (R/z)) 



8i:r]R 3 



(14) 



with r] denoting the viscosity of the fluid. (3 rr (l/z) and 
P dr (R/z) are dimensionless functions including the z- 
dependence. A numerical scheme that allows to accurately 
calculate these functions at arbitrary heights can be found 
in [16, 17]. In order to compare the results of the model 
presented here with the numerical solutions of the MFPT 
problem we use mean-field values in the following way 

oo oo 

Dg:= [ dz^ s (z)Dg{z), Ag:= f dz* s (z)Ag(z). (15) 



That is we take for Dg and Ag averages with respect to the 
stationary probability function ^ s (z) from eq. (3) of the 
z-dependent quantities Dg(z), Ag(z) defined in eq. (14). 

Comparison to simulation results. In fig. 2 wc 
compare our results for the angle averaged MFPT (T)g> 
as obtained from the analytical calculation (solid lines) 
and from computer simulations (symbols). For this pur- 
pose we rescale time in units of the diffusive time scale 
6TTTIR 3 /ksTa. In fig. 2a (T)$i is shown as a function of 
the Pcclct number Pe := 6nriR 3A //kBT a (other parame- 
ters are defined in the figure caption) . The Peclet number 
is a dimensionless measure for the relative importance of 
deterministic and diffusive motion and Ag oc Pe. In the 
case of zero rotational drift, Ag = 0, Pe = and motion 
is purely diffusive. In addition to the simulation results 
(symbols) fig. 2a also shows our main result (lines) eq. (12) 
where we used eqs. (13) and (15) for the rate k and the 
constants Ag,Dg, respectively. The agreement between 
the mean-field and simulation results is surprisingly good 
for the parameter values chosen. In particular, even the 
small shoulder for N r — 1 seems to be reproduced by the 
analytical result. In fig. 2b (T)g> is shown as a function of 
the number of receptor patches N r in the diffusive limit 
Ag = and for two different values of the patch height ro . 
Whereas the theoretical approximation with the rate given 
by eq. (13) works quite well for r = 10~ 3 R, for a tenfold 
larger ro we notice clear deviations. The larger ro the 
larger 9q and the assumption of the dwelling time within 
[0, 26>o] being small is less valid. Nevertheless, eq. (12) de- 
scribes the functional dependence of {T(N r ))g> in a quali- 
tative way. In fact quantitative agreement can be achieved 
by using the reaction rate A: as a fit parameter (not shown). 

Asymptotic limits. — The results shown in fig. 2 
were obtained for a fixed vertical drift with strength given 
by Pe z — 50. Here, Pe z :— F z R/k B T a is the Peclet num- 
ber for the z-direction. Fig. 3a displays the comparison be- 
tween the numerical result of {T)g> and the approximation 
eq. (12) as a function of Pe z in the diffusive limit Ag w 0. 
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Fig. 2: Comparison between numerical results for the MFPT 
(T)gi averaged over all initial orientations and the mean- field 
approximation according to eqs. (12,13,15). The numerical re- 
sults where obtained as described in ref. [13] for an initial height 
2 = 2R. In order to display the MFPT with initial height 
z — R + ro the falling time T(R + ro\2R) was subtracted ac- 
cording to eq. (1). (a) (T)gi as a function of the dimensionless 
Peclet number Pe for different numbers of receptor patches 
N r . (b) (T)gi as a function of N r for two different values of the 
patch height ro. (Other parameters: patch radius r v = 0.001.R, 
vertical drift F z = Pe z R/k B T a with Pe z = 50.) 



Fig. 3: (a) The numerically obtained MFPT (T)g> as a function 
of the strength of the vertical drift Pe z (+). The dashed lines 
show the asymptotic behaviour (T)gi oc 1/Pe z for small Pe z 
and the value of {T)g> around Pe z = 10 3 . The solid line shows 
the mean-field approximation according to eqs. (12,13,15). 
(N T = 10, r = r p = 10' 3 R). (b) The different scaling be- 
haviour of (T)gi in regard to the number of receptor patches 
N r for zero and large rotational drift Ag at Pe z = 10 3 is shown. 
The theoretical result from eq. (12) matches well the numerical 
results (symbols) when for k infinity and k = 2Ag is chosen for 
Pe = and Pe « 1500, respectively (r = r p = 10~ 3 R). 



For small Pe z the numerical results for the MFPT shows 
a 1/ Pe z scaling behaviour as indicated by the dashed line 
in fig. 3a. For larger Pe z , (T)gi plateaus. At even larger 
Pe z the MFPT increases with increasing Pe z (not shown) 
as Dg(z) — > for z — > R [17]. One sees that the approx- 
imation eq. (12) underestimates the numerical result for 
small values of Pe z and overestimates it for large Pe z . For 
intermediate values of Pe z of the order of 10-100 we find 
good agreement between the two results as demonstrated 
before in fig. 2. For large Pe z w 10 3 the numerical result 
provides (T)g> w 0.1. About the same value is given by 
the first term in eq. (12). This means that in the diffu- 
sive limit and for large Pe z an encounter occurs almost 
with probability one for 9 £ [0, 20n]- This is plausible as 
then the duration time for 9 £ [0, 2f?n] is long enough such 
that the particle will most probably encounter a height 



z < R + ro while the receptor patch points downwards. 
Therefore, the rate k for Pe z « 10 3 is rather infinity than 
the value given by the estimate eq. (13). At small values of 
Pe z we expect combining eq. (12) and eq. (13) the second 
term in eq. (12) to be dominant and (T)$i oc l/(DgPe z ). 
As the diffusion coefficient Dg(z) is a monotonically in- 
creasing with increasing distance z from the wall also the 
term Dg defined in eq. (15) becomes larger with decreas- 
ing Pe z . This is the reason why the estimate eq. (12) for 
(T)gi does not provide the right scaling behaviour. On 
the other hand using De|pe z =50 also at larger Pe z values 
a much better agreement (not shown) between the theo- 
retical estimate and the numerical result of the MFPT is 
obtained. That implies that the faster rotations far away 
from the wall, which are included in the definition eq. (15), 
are not relevant for the encounter process which happens 
only very close to the wall. 
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In fig. 3b we demonstrate that eq. (12) matches well 
to the numerically obtained MFPT for Ag = and large 
Pe z when k is set to infinity as discussed above. There, 
for Pe z = 10 3 , the numerical result for (T)g> (+) and the 
theoretical result eq. (12) (full line) are shown. For the 
dependence on the number of receptor patches N r we hnd 
in this limit (and for N r small) 

A9 3 4tt 2 

{T)e ' ~ m^e * T^W' (16) 

i.e. the MFPT approximately scales as The situ- 

ation is different in the deterministic limit, i.e. for large 
Ag. Then, eq. (12) is approximately 




i.e. the MFPT scales as 1/N r at small numbers of receptor 
patches N r . In fig. 3b we also show a comparison between 
eq. (17) (dashed line) and the numerically obtained MFPT 
(x). Here we find that good agreement is obtained when 
we choose k « 2Ag, which is of the same order of magni- 
tude as the estimate given by eq. (13). Thus, in contrast to 
the diffusive limit, the second and the first term in eq. (12) 
are of the same order. 

Summary and outlook. — In this letter we have 
derived an approximate expression eq. (12) for the angle 
averaged MFPT for receptor-ligand encounter between a 
sphere cquidistantly covered with receptor patches and a 
wall homogeneously covered with ligands in a 2D geom- 
etry. The main idea of our analysis was to integrate out 
the motion in z-direction by absorbing falling in a reaction 
term for the rotational degree of freedom. The coefficients 
for the diffusion and drift terms of this equation were esti- 
mated from mean field arguments in eq. (15). Our deriva- 
tion was based on the central assumption that the heights 
of the sphere at two successive times of a receptor patch 
pointing downwards is uncorrelated. This is not true in 
general, and accordingly the result derived for the reaction 
rate k in eq. (13) is valid only for a small range of param- 
eter values. However, the rate k can also be viewed as a 
fit parameter. In that case eq. (12) matches the results 
obtained in computer simulations of the full problem over 
a large range of parameters. 

In the future, our analysis could be extended in different 
ways. Here we have only considered a homogeneous lig- 
and density. For non-homogeneous ligand coverage pi < 1 
not every receptor- wall encounter is a productive receptor- 
ligand encounter. This might again be expressed by an 
appropriate choice of the rate k with k oc p\. Even for 
homogeneous ligand density, not every receptor-ligand en- 
counter has to lead to functional adhesion under flow. 
Depending on the receptor-ligand system under considera- 
tion, additional steps might be required to achieve a stable 
bond. Conceptually, one could regard the first encounter 
as formation of an encounter complex [2]. The next step 



would then be the transition into a final complex. In a 
such a two-state system, also unbinding becomes impor- 
tant, as it characterizes the stability of the final complex. 
Again this process might be expressed by an appropriate 
choice of the rate k as used here. If complete unbinding 
occurs, there is also the possibility that the cell forms a 
new bond downstream of the old one. At sufficient den- 
sities of receptors and ligands, this mechanism eventually 
leads to the physiologically very relevant process of rolling 
adhesion. Recently we have extended the computer sim- 
ulations used above to measure MFPTs to also simulate 
the process of rolling adhesion [18]. Future work has to 
show how these simulations can now be made more effi- 
cient using the dimensional reduction introduced here. In 
general, analytical progress by dimensional reduction and 
introduction of appropriate reaction terms might be a very 
promising strategy also for other biological systems which 
involve a complex interplay between different transport 
modes. 
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